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Table 1 (conn 

...readjmps 

if 0 < 0 && dxfdmax] .offset &&fj){ 

(void) Iseek(fd, dx[dmax] .offset, 0); 
(void) read(fd, (char *)&dx[dmax].jp, sizeof (struct jmp)); 
(void) read(fd, (char *)&dx[dmax].offset, sizeof(dx[dmax]. offset)); 
dx[dmax].ijmp = MAXJMP-1; 

> 

else 

break; 

} 

if (i > = JMPS) { 

fprintf(stderr, " %s: too many gaps in alignment\n", prog); 
cleanup(l); 

} 

if(j > = 0){ 

siz = dx[dmax].jp.nD]; 
xx = dx[dmax].jp.x|j]; 
dmax + = siz; 

if (siz < 0) { /* gap in second seq */ 

pp[l].n[il] = -siz; 
xx += siz; 

/* id = xx - yy + lenl - 1 
*/ 

pp[l].x[il] = xx - dmax + lenl - 1; 
gapy++; 

i ngapy -= siz; 

J H /* ignore MAXGAP when doing endgaps */ 

7 siz = (-siz < MAXGAP 1 1 endgaps)? -siz : MAXGAP; 

p30 > 

U else if (siz > 0) { /* gap in first seq */ 

m PP [0].n[iO] = siz; 

?i pp[0].x[i0] = xx; 

uU s gapx++; 
r ^ ngapx + = siz; 

/* ignore MAXGAP when doing endgaps */ 

siz = (siz < MAXGAP 1 1 endgaps)? siz : MAXGAP; 

i0++; 

40 > } 

else 

break; 

> 

45 /* reverse the order of jmps 

*/ 

for 0 = 0, i0~; j < iO; j++, iO-) { 

i = pp[0].nQ]; pp[0].n[j] = pp[0].n[i0]; pp[0].nTi0] = i; 
5Q i = PP[0].x[j]; pp[0].x[j] = pp[0].x[i0]; pp[0].x[i0] = i; 

for(j =0, il-;j < il;j + +,il-){ 

i = pp{l].nO]; pp[l].nrj] = pp[l].n[il]; pp[l].n[il] = i; 
^ i = pp[l].xD3; pp[l].xQ] = pp[l].x[il]; pp[l].x[il] = i; 

55 if (fd > = 0) 

(void) close(fd); 

if(ffi{ 

(void) unlink(jname); 

60 offset '= 0; 

} 

} Page 3 of nwsubr.c 



151 



Table 1 (conV) 

/* 

* write a Med jmp struct offset of the prev one (if any): nw() 
*/ 

5 writejmps(ix) writeimps 
int ix; 

{ 

char *mktempO; 
10 if(!fj){ 

if (mktempQname) < 0) { 

fprintf(stderr, "%s: can't mktempO %s\n", prog, jname); 
cleanup(l); 

1« } 

1 -> if ((fj = fopen(jname, "w")) = = 0) { 

fprintf(stderr, "%s: can't write %s\n\ prog, jname); 
exit(l); 

} 

} 

20 (void) fwrite((char *)&dx[ix].jp, sdzeof(struct jmp), 1, fj); 

(void) fwrite((char *)&dx[ix] .offset, sizeof(dx[ix].offset), 1, fj); 
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Table 2 



PRO XXXXXXXXXXXXXXX (Length = 15 amino acids) 

Comparison Protein XXXXXYYYYYYY (Length = 12 amino acids) 

5 % amino acid sequence identity = 

(the number of identically matching amino acid residues between the two polypeptide sequences as determined 
by ALIGN-2) divided by (the total number of amino acid residues of the PRO polypeptide) = 

10 5 divided by 15 = 33.3% 
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